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Summary. — We present an experimental study of the time of flight properties 
of a gas of ultra-cold fermions in the BEC-BCS crossover. Since interactions can 
be tuned by changing the value of the magnetic field, we are able to probe both 
non interacting and strongly interacting behaviors. These measurements allow us 
to characterize the momentum distribution of the system as well as its equation 
of state. We also demonstrate the breakdown of superfluid hydrodynamics in the 
weakly attractive region of the phase diagram, probably caused by pair breaking 
taking place during the expansion. 



1. — Introduction 

Feshbach resonances in ultra cold atomic gases offer the unique possibility of tun- 
ing interactions between particles, thus allowing one to study both strongly and weakly 
interacting many-body systems with the same experimental apparatus. A recent major 
achievement was the experimental exploration of the BEC-BCS crossover [TJ [2j [3l HJ El E] , 
a scenario proposed initially by Eagles, Leggett, Nozieres and Schmitt-Rink to bridge the 
gap between the Bardeen-Cooper-Schrieffer (BCS) mechanism for superconductivity in 
metals, and the Bose-Einstein condensation of strongly bound pairs [TJ [9] . Here, we 
present a study of the crossover using time of flight measurements. This technique gives 
access to a wide range of physical properties of the system and has been successfully 
used in different fields of physics. The observation of elliptic flows was for instance used 
to demonstrate the existence of quark-gluon plasmas in heavy ion collisions [TP] , In cold 
atoms, the ellipticity inversion after free flight is a signature of Bose-Einstein condensa- 
tion [TTJ [12] . In an optical lattice the occurrence of interference peaks can be used as the 
signature of the superfluid to insulator transition [13] and, with fermions, it can be used 
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to image the Fermi surface [14 . Two series of time-of- flight measurements are presented: 
expansion of the gas without interactions, which gives access to the momentum distri- 
bution, a fundamental element in the BEC-BCS crossover, or with interactions, which 
allows us to characterize the equation of state of the system, and probe the validity of 
superfluid hydrodynamics. 

2. — Experimental method 

In a magnetic trap, a spin polarized gas of N = 10 6 6 Li atoms in \F = 3/2, mp = 
+3/2) is sympathetically cooled by collisions with 7 Li in \F = 2, tuf = +2) to a tempera- 
ture of 10 /iK. This corresponds to a degeneracy of T/Tp ~ 1, where T F = huu^N) 1 / 3 /k B 
is the Fermi temperature of the gas. The magnetic trap frequencies are 4.3 kHz (76 Hz) 
in the radial (axial) direction, and Q = {uJxUJyUz) 1 ^ is the mean frequency of the trap. 
Since there are no thermalizing collisions between the atoms in a polarized Fermi gas, the 
transfer into our final crossed dipole trap, which has a very different geometry (FlG. [T]), 
is done in two steps. We first transfer the atoms into a mode-matched horizontal sin- 
gle beam Yb:YAG dipole trap, with a waist of ~ 23 /im. At maximum optical power 
(2.8 W), the trap depth is ~ 143 /iK (15 Tp), and the trap oscillation frequencies are 
6.2(1) kHz (63(1) Hz) in the radial (axial) direction. The atoms are transferred in their 
absolute ground state \F = 1/2, mp = +1/2) by an RF pulse. We then sweep the mag- 
netic field to 273 G and drive a Zeeman transition between \F — 1/2, mp = +1/2) and 
\F = 1/2, ttif — —1/2) to prepare a balanced mixture of the two states (better than 5%). 
At this magnetic field, the scattering length between both states is —280 ao (Fig. [2j). Af- 
ter 100 ms the mixture has lost its coherence, initiating collisions in the gas. During the 
thermalization process about half of the atoms are lost. We then perform a final evapo- 
rative cooling stage by lowering the trap depth to ~ 36 /iK. At this point, we ramp up 
a vertical Nd:YAG laser beam (power 126 mW and waist ~ 25 /im), obtaining our final 
crossed dipole trap configuration (Fig. [1]). The measured degeneracy is T/Tp < 0.15. 
The magnetic field is then increased to 828 G (in the vicinity of the Feshbach resonance, 
see Fig. [2]), where we let the gas thermalize for 200 ms before performing subsequent 
experiments. 

3. — Momentum distribution 

In standard BCS theory, the ground state of an homogeneous system is described by 
a pair condensate characterized by the many-body wave function 

k 

where |0) is the vacuum and a\, is the creation operator of a fermion with momentum 
k and spin a. In this expression, \vk\ 2 can be interpreted as the occupation probability 
in momentum space, and is displayed in Fig. [3^l for several values of the interaction 
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Fig. 1. - Ioffe-Pritchard trap and crossed 
dipole trap used for the experiments. The 
crossed geometry allows us to change the as- 
pect ratio of the trap. 




Magnetic Field [G] 

Fig. 2. - 6 Li Feshbach resonance between \F — 
1/2, m F = +1/2) and \F = 1/2, m F = -1/2). 
The broad Feshbach resonance is located at 
834 G. The balanced mixture is prepared at 
273 G. 



parameter I/Zcf^, where k F is the Fermi wave vector of the non interacting Fermi gas 
(E F = H 2 k 2 F /2m). One effect of the pairing of the atoms is to broaden the momentum 
distribution. In the BCS limit (1/fc^a — > — oo), the broadening with respect to the 
momentum distribution of an ideal Fermi gas is very small, of the order of the inverse of 
the coherence length £. In the unitary limit (l/kpa — > 0) it is expected to be of the order 
of kp. In the BEC limit {l/kpa — > oo) we have molecules of size a so the momentum 
distribution, which is given by the Fourier transform of the molecular wave function, has 
a width h/a. 

In a first series of expansion experiments, we have measured the momentum distribu- 
tion of a trapped Fermi gas in the BEC-BCS crossover. Similar experiments have been 
performed at JILA on 40 K [T5] . 

In order to measure the momentum distribution of the atoms, the gas must expand 
freely, without any interatomic interactions. To achieve this, we quickly switch off the 
magnetic field so that the scattering length is brought to zero (see Fig. [2]) [16]. We 
prepare N = 3 x 10 4 atoms at 828 G in the crossed dipole trap with frequencies uo x = 
2tt x 2.78 kHz, u y = 2tt x 1.23 kHz and u z = 2tt x 3.04 kHz. The magnetic field 
is adiabatically swept in 50 ms to different values in the crossover region. Then, we 
simultaneously switch off both dipole trap beams and the magnetic field (with a linear 
ramp of 296 G//is). After 1 ms of free expansion, the atoms are detected by absorption 
imaging. The measured density profiles give directly the momentum distribution of the 
gas integrated along the imaging direction. 

In Fig. [3j we show the measured momentum distributions for three different interac- 
tion parameters, corresponding to the BCS side of the resonance, the unitary limit and 
the BEC side of the resonance. Together with our data, we have plotted the predictions 
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of mean field BCS theory at T = 0, taking into account the trapping potential with 
a local density approximation [17J. k° F is now the Fermi wave- vector calculated at the 
center of the harmonic trap for an ideal gas. 

Some precautions need to be taken concerning this type of measurements due to 
possible density dependent losses during the magnetic field switch-off. If the magnetic 
field is not turned off fast enough, some atoms can be bound into molecules while the 
Feshbach resonance is crossed. The molecules are not detected with the imaging light 
and therefore will appear as a loss of the total number of atoms. Even if, as in our case, 
the Feshbach resonance is crossed in only 1 /is, this time may not be small compared to 
the typical many-body timescale (h/Ep ~ 1.3 jis for Fig. 3 data). 

To investigate quantitatively this effect, we have performed an additional experiment 
in a more tightly confining trap. We prepare a gas of 5.9 x 10 4 atoms at 828 G in a trap 
with frequencies u x = 2tt x 1.9 kHz, uo y = 2tt x 3.6 kHz and uj z = 2tt x 4.1 kHz. The total 
peak density in the trap is 1.3 x 10 14 atoms/cm 3 . We let the gas expand at high field 
for a variable time ts, then switch off B and detect the atoms after 0.5 ms of additional 
free expansion. Assuming hydrodynamic expansion at unitarity we calculate the density 
after ts [18J and obtain the fraction of atoms detected as a function of the density of 
the gas when the resonance is crossed. For instance, we detect c± 60% fewer atoms for 
tg = compared to ts = 0.5 ms, where the density is a factor ~ 10 3 lower. The results 
are nicely fitted by a Landau- Zener model : 

^detected/^total = exp f-A^^ 



2B 



where n(t#) is the total density at B the sweep rate and A the coupling constant 
between the atoms and the molecules. We determine A ~ 5 x 10 -12 G m 3 /s. Our result 
is five times smaller than the MIT value A ~ 24 x 10 -12 G m 3 /s [19], measured at a total 
peak density of 10 13 atoms/cm 3 (one order of magnitude smaller than in our experiment). 
The theoretical prediction, assuming only two body collisions, is A = 19 x 10 -12 G m 3 /s 
[20] . The difference between our measurement and theory may suggest that many-body 
effects are important in our case. Finally, using our value of A the model predicts an 
atom number loss of about 27% for the momentum distribution measurements of Fig. 3. 
This loss is comparable to our shot-to-shot fluctuations in atom number and therefore 
was not unambiguously observed. 

In conclusion, we have performed a measurement of the momentum distribution of 
a trapped Fermi gas. The results are found in reasonable agreement with BCS theory 
despite the fact that it is not expected to be quantitatively correct in the strongly in- 
teracting regime. In future work, experiments at lower density will be performed, in 
order to avoid the observed loss effect. This should allow us to distinguish between BCS 
and more exact theories [2T] , It would also be interesting to perform measurements at 
different temperatures as in Ref. [22] . 
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Fig. 3. - (a): Momentum distribution of a uniform Fermi gas for 1/kFd = —1 (solid line), 
1/kFCL = (dotted line) and 1/kFd = 1 (dashed line) calculated from mean field BCS theory 
at T — [17]. The results obtained from quantum Monte Carlo simulations [21 show that 
BCS theory slightly underestimates the broadening; (b): Measured momentum distribution of 
a trapped Fermi gas on the BCS side of the resonance (l/k F a = —0.42); (c): Unitary limit 
(l/k F a = 0); (d): BEC side of the resonance (l/k F a — 0.38). The solid lines in (b), (c) and (d) 
are the predictions of BCS mean field theory taking into account the trapping potential with a 
local density approximation [IT]. k F is defined in the text. 



4. — Release energy 

In a second series of experiments, we have performed expansions at constant magnetic 
field, thus keeping the interactions present during the time of flight. The analysis of size 
measurement across the BEC-BCS crossover yields valuable information on the influence 
of interactions on the properties of the system. In particular, we have measured the 
release energy of the gas in the BEC-BCS crossover [3 . On resonance (1/k^a = 0), the 
gas reaches a universal behavior [23] . The chemical potential fi is proportional to the 
Fermi energy \i = (1 + f3)Ep. We have determined the universal scaling parameter f3 
from our release energy measurement. 
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The starting point for the experiment is a nearly pure molecular condensate of 7 x 10 4 
atoms at 770 G, in an optical trap with frequencies uj x = 2tt x 830 Hz, uj y = 2ty x 2.4 kHz, 
and ou z = 2tt x 2.5 kHz. We slowly sweep the magnetic field at a rate of 2 G/ms to various 
values across the Feshbach resonance. We detect the integrated density profile after a 
time of flight expansion of 1.4 ms in several stages: 1 ms of expansion at high magnetic 
field, followed by a fast ramp of 100 G in 50 [is in order to dissociate the molecules and, 
after the fast switch-off of the magnetic field, 350 /as of ballistic expansion. 

Fig. [4] presents the gas energy released after expansion, which is calculated from 
gaussian fits to the optical density after time of flight: E' rel = m(2a y + <t^)/2t 2 , where 
<Ji is the rms width along i and r is the time of flight. We assume that the size a z 
(which is not observed) is equal to a y . Note that both in the weakly interacting case 
and unitarity limit the density has a Thomas-Fermi profile and the release energy can 
be calculated from the exact profiles. However, we have chosen this gaussian shape to 
describe the whole crossover region with a single fit function. This leads to a rescaling 
of the release energy. In particular, the ideal Fermi gas release energy in an harmonic 
trap is E re \ = 3/SEp but when using the gaussian fit to the Thomas- Fermi profile we get 
instead E' TeX = 0A6E F as shown in Fig. HJ 

The release energy in the BEC-BCS crossover varies smoothly. It presents a plateau 
for —1/kpa < —0.5, (BEC side) and then increases monotonically towards that of 
a weakly interacting Fermi gas. On resonance, the release energy scales as E re \ = 
yjl + (3 E® el , where E® el is the release energy of the non interacting Fermi gas. The 
square root comes from the average over the trap. At 834 G, we get (3 = —0.59(15). 
This value is slightly different from our previous determination (3 = —0.64(15), where 
the resonance was assumed to be located at 820 G instead of 834 G [3]. Our result 
agrees with other measurements performed on 6 Li and with some theoretical predictions 
(see Table [J). Remarkably, the recent 40 K measurement at JILA is also in very good 
agreement, thus proving the universality of the unitarity regime. 



5. - Ellipticity 

Nontrivial information can be extracted from the measurement of the aspect ratio 
of the cloud after expansion. For instance, in the first days of gaseous Bose-Einstein 
condensates, the onset of condensation was characterized by an ellipticity inversion after 
time of flight, a dramatic effect compared to the isotropic expansion of a non condensed 
Boltzmann gas. In the case of strongly interacting Fermi gases, ellipticity measurements 
can be used as probes for the hydrodynamic behavior of the system, and constitute an 
indirect signature of the appearance (or breakdown) of superfluidity. 

We have studied the ellipticity of the cloud as a function of the magnetic field for 
different temperatures. As before, the density profiles are fitted with gaussians, and the 
ellipticity is defined as r] = cr y /a x . We prepare iV = 3x 10 4 atoms at 828 G in a crossed 
dipole trap. The magnetic field is adiabatically swept in 50 ms to different values in the 
crossover region. Then, we switch off both dipole trap beams and let the gas expand 
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Fig. 4. - Rescaled release energy E' Tel of a trapped Fermi gas in the BEC-BCS crossover as a 
function of the magnetic field and as a function of —l/k° F a [3]. The dashed line is the rescaled 
release energy of a T — non interacting Fermi gas. From the measurement at resonance we 
extract = -0.59(15). 
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Experimental results on 6 Li 
at finite T 


This work 
ENS 2004 [3] 
Innsbruck [24] 
Duke [25] 
Rice [26] 


-0.59(15) 
-0.64(15) 

— D 7^+ 12 
u - ' °-0.09 

-0.49(4) 
-0.54(5) 


Experimental result on 40 K 
extrapolation to T=0 


JILA [27] 


-0.54lg;E 


Theoretical predictions 
at T = 


BCS theory [3 E] 
Astrakharchik et al. [28] 
Carlson et al. [29] [30] 
Perali et al. [31] 
Pade approximation [23 , 32 
Steel [23] 

Haussmann et al. [3l] 


-0.41 

-0.58(1) 

-0.58(1) 

-0.545 

-0.67 

-0.56 

-0.64 


Theoretical predictions 
at T — T c 


Bulgac et al. [35] 
Burovski et al. [36] 


-0.55 
-0.507(14) 



Table I. - List of the recent experimental measurements and theoretical predictions of the uni- 
versal scaling parameter p. 
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for 0.5 ms in the presence of the magnetic field. After 0.5 ms of additional expansion 
at B = 0, the atoms are detected by absorption imaging. Fig. [5^ and Fig. [5b show 
the measured value of the ellipticity as a function of the magnetic field for two different 
samples, which are at different temperatures. Together with the experimental results we 
have plotted the expected anisotropy from superfluid hydrodynamics [18]. For this, we 
have extracted from the quantum Monte Carlo simulation of ref. [28J the value of the 
poly tropic exponent 7, defined as 7 = ~f^- 

The first series of measurements is done in a trap with frequencies uo x = 2tt x 1.39 kHz, 
ujy = 2tt x 3.09 kHz, co z = 2tt x 3.38 kHz and trap depth ~ 1.8 Tp. The measured ellipticity 
(Fig. [5^) is in good agreement with the hydrodynamic prediction on the BEC side, at 
resonance and on the BCS side until l/k^a = —0.15. It then decreases monotonically to 
1.1 at l/k° F a = -0.5. 

For the second series of experiments we prepare a colder sample in a trap with fre- 
quencies uj x = 2tt x 1.24 kHz, uo y = 2tt x 2.76 kHz, uj z = 2ty x 3.03 kHz and trap depth 
~ 1.6 Tp. In this case the behavior of the anisotropy is very different (Fig. [5}}). We 
observe a plateau until 1/k^a = —0.33, in good agreement with the hydrodynamic pre- 
diction, and at this critical magnetic field there is a sharp decrease of 77 to a value close 
to 1.2. This sharp transition seems analogous to the sudden increase of the damping of 
the breathing mode observed in Innsbruck [37] . 

In a third experiment, we measure the ellipticity at unitarity as a function of trap 
depth (hence of the gas temperature). Below a critical trapping laser intensity, the 
ellipticity jumps from a low value (1.1) to the hydrodynamic prediction 1.45. 

In all cases, the decrease of the anisotropy indicates a breakdown of superfluid hydro- 
dynamics in the weakly attractive part of the phase diagram or at higher temperature. 
A first possibility would be that the gas crosses the critical temperature in the trap. 
However, we know from the MIT experiment [38] that pair breaking can occur during 
the expansion. During the time of flight, both the density and kp decrease. On the BEC 
side of resonance, the binding energy of the molecules (—h 2 /ma 2 ) does not depend on 
the density and the pairs are very robust. By contrast, on the BCS side of resonance 
the generalized Cooper pairs become fragile as the gap decreases with 1/kpa and they 
can be broken during the expansion. Our experiments use the ellipticity of the cloud as 
a probe and are complementary to the MIT approach, where the breakdown of superflu- 
idity was characterized by the disappearance of vortices during the expansion of the gas. 
We are planning additional experiments in order to investigate wether the breakdown of 
superfluidity occurs in the trap or during the expansion. 



6. — Conclusion 

The results presented here constitute a first step in the understanding of the free flight 
properties of strongly correlated fermionic systems. In future work, we will investigate 
more thoroughly the pair breaking mechanism taking place during the expansion in the 
BCS part of the phase diagram. We point out the need for a dynamic model of the 
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Fig. 5. - Ellipticity of the gas after expansion from a trap of depth ~ 1.8 Tf (a) and from a trap 
of depth rsj 1.6 Tf (b). Solid lines: hydrodynamic predictions. 



expanding gas at finite temperature. 
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